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Abstract 

We present an open-loop (bang-bang) scheme which drives an open two-level quantum system 
to any target state, while maintaining quantum coherence throughout the process. The control is 
illustrated by a realistic simulation for both adiabatic and thermal decoherence. In the thermal 
decoherence regime, the control achieved by the proposed scheme is qualitatively similar, at the 
ensemble level, to the control realized by the quantum feedback scheme of Wang, Wiseman, and 
Milburn [Phys. Rev. A 64 063810 (2001)] for the spontaneous emission of a two-level atom. The 
performance of the open-loop scheme compares favorably against the quantum feedback scheme 
with respect to robustness, target fidelity and transition times. 
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I. INTRODUCTION 



Over the last few years, a number of open-loop (bang-bang) control schemes have been 
proposed to eliminate the effects of decoherence for a single two-level quantum system in 
contact with the environment, by using an external control DDQD . This type of control, 
pioneered by Lloyd and Viola (see Refs. l], SQ] arid references therein), relies on applying 
frequent control pulses to the two-level system in order to cancel the decoherence effects of 
the system-environment interaction. Vitali and Tombesi ?4| have also considered applying 
a sequence of frequent parity kicks, as well as an appropriate stochastic modulation, to 
achieve the same goal. It has been shown that, in either approach, the decoherence of the 
system can be effectively suppressed if the pulse rate is much higher than the decoherence 
rate due to the system-environment interaction. In other words, a fast rate control "freezes" 
decoherence in a manner similar to the quantum Zeno effect. 

Alternative approaches have also been considered to eliminate or mitigate the undesirable 
effects of decoherence in open quantum systems, e.g., decoherence free subspaces (DFS) 
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quantum error correction and quantum feedback 

Inspired by Lloyd and Viola's seminal results, we have recently proposed a new imple- 
mentation of open-loop quantum control [12], which tailors control pulses more efficiently, 
by taking into account the possible knowledge of the decoherence function. Several new 
features of the scheme are: (i) decoherence and control act simultaneously within a realistic 
model that includes them from first principles; (ii) the required control is directly related 
to and calculated from the decoherence effects, which presents the practical advantage of 
maintaining the frequency and amplitude of the required controls at minimal levels; and 
(iii) the effect of imperfect pulses on the efficiency of the control was assessed and found 
to be quite tolerable, even for rather large noise amplitudes. In this paper, our aim is to 
minimize the effect of decoherence due to interactions between the quantum system and its 
environment, while concommitantly driving the system from an initial state to an arbitrary 
target state. 

Recently, KKeeeJa et al. Q used M contro, to steer a qU aat Um SyS te m to a ta rg et 
state in a minimum time, while assuming that the effects of relaxation are negligible over 
that timescale. They have also demonstrated that controllability for a closed quantum 
system is equivalent to the universality that is demanded from a quantum computer. In 
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follow-on work, Khaneja et al. [14] used optimal control to minimize the difference between 
the actual density matrix at the final time and the target density matrix, for two spins in 
the presence of relaxation. We note that their relaxation model relies on phenomenological 
decay rates to describe general decoherence processes, without any explicit reference to the 
environment, and the control is applied a posteriori, to the "reduced" equations. This is 
in contrast with the present scheme, where the interaction with the environment and the 
control are accounted for from first principles, in the full Hamiltonian description of the 
system. 

We illustrate our approach for pure states of a two-level system, i.e., a qubit, which is 
in contact with the environment. Both the initial and the target states are known a priori, 
and the unitary evolution of the two-level system is driven by frequent control pulses that 
are customized, if the decoherence function is known. This is the case - at least in principle 
- if one has an explicit model for the environment and its interaction with the qubit. 

Using numerical simulations, we show explicitly that the qubit can be driven between any 
two points (initial and target states) on the curve determined by the intersection of either 
the x — z, or the y — z planes with the Bloch sphere. In general, the state of the qubit can 
be driven between any two points on the surface of the Bloch sphere by using a sequence of 
two different control Hamiltonians. 

In one particular instance, namely for thermal decoherence and external control propor- 
tional to the o y polarization, the system we consider is formally similar to the system used 
in the quantum feedback scheme proposed by Wang et al. [10], which models the evolu- 
tion of a two-level atom that emits spontaneously, using stochastic quantum trajectories. 
Photocurrent feedback from perfect homodyne detection of spontaneous emission is used to 
alter the atom-environment dynamics, to drive the atom to a pure target state. Of course, 
the feedback delay time has to be shorter than the spontaneous emission timescale, in order 
to successfully maintain the atom in the target state. 

Whenever warranted, we shall discuss the performance of the proposed open-loop control 
in comparison with the performance of quantum feedback control, in order to understand 
the limitations of each scheme, and characterize the regimes in which they work best. In so 
doing, we also generalize the original feedback control scheme propsed by Wang et al., by 
using different control Hamiltonians. 

The remainder of the paper is organized as follows: In Section 2, we give an overview of 
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the spin-boson model used in our targeting control scheme, and in Section 3 we outline the 
control strategy employed to achieve a particular target state. In Section 4, we briefly review 
the quantum feedback scheme proposed by Wang et al. |l0| for a spontaneously-emitting 
two-level atom. Finally, in Section 5, we present the results of our numerical simulations, 
and discuss them in comparison with the results of the quantum feedback scheme. 



II. THE MODEL 

The spin-boson model |l5] is often used to describe a two- level system {|1 >, |2 >} 
interacting with a large number of boson modes representing the environment. The complete 
Hamiltonian 

H = H S + H e + H l + H c , (1) 

accounts for the two- level system, the environment, the interaction between the two, and an 
external control. The Hamiltonian of the two-level system is 

2 

H s = J2 p u E i, (2) 

i=i 

where Pa denote the projection operators Pa = and Ei are the corresponding energies. 

We take E\ = —Kujq/2 and Ei = huo/2, so that H s = hujQa z . The Hamiltonian of the 
environment is 

CO 

H e = h^2uJo g ala q , (3) 

q=l 

where a q , a) are the annihilation and creation operators respectively, for the g-th boson 
mode, and uiQ q is the corresponding frequency. 

The interaction between the two- level system and the environment, 

oo 

Hi = -he(a x a x + f3 z a z ) ^(fl*a 5 + fi g aj) (4) 

q=l 

is responsible for the quantum decoherence. Adiabatic decoherence (a x — 0, /3 Z — 1) acts 
on a short timescale, leading to the decay of the off-diagonal density matrix elements. On 
a longer timescale, thermal decoherence (a x — 1, f3 z — 0) changes all of the density matrix 
elements, and leads to the exponential decrease of the excited state population of the two- 
level system. These two timescales can be accounted for by adjusting the relative magnitude 
of the coefficients, a x and (3 Z . The interaction with the environment is parametrized by e 



(whose magnitude indicates the strength of the coupling) and can result in a phase flip or 
a bit flip (a bit-phase flip can also be obtained either by combining the two effects or by 
including a a y term). However, from a formal viewpoint, the effects of a x and a y in the 
interaction Hamiltonian are essentially similar, and we shall not consider the latter as far 
as Hi is concerned. 

In general, the control Hamiltonian is proportional to a linear superposition of the pseudo- 
spin operators a x and a y , 

H c = -hQ F V(t)[c x a x + CyO-y], (5) 

which can be specified by adjusting the phase of the driving field. The control is a classical 
coherent driving field, which is applied as a series of frequent pulses, V(t). We assume that 
these pulses act strongly upon the two-level system, via its polarization, but do not affect 
the environment significantly. 

The dynamics of the two-level system is expressed in terms of parameters normalized to 
the Rabi frequency, Qp. The time evolution is described in terms of dimensionless Rabi time 
units, r = Q F t, and all the frequencies, namely u q = u oq /Qp, cu = u /Qp, g q = fl q /flp and 
U12 = CIu/Qf, are also renormalized with respect to the Rabi frequency. The frequency Qu 
is the transition frequency between the energy levels of the two-level system. 

To simplify the calculations and render them more transparent, it is convenient to write 
the evolution of the system in the interaction representation, and consider the adiabatic 
and thermal decoherence regimes separately. The unitary evolution operator in the inter- 
action representation, Ui, is related to the unitary evolution operator in the Schrodinger 
representation, U, by 

U I (s, e, r) = e iHoT/n U(s, e, r) (6) 

where H = H s + H e . 

In this representation, the evolution of the density matrix of the whole system (qubit 
plus environment) reads 

where the interaction Hamiltonian contains both the interaction with the environment and 
the applied control: 

Hj = H Ic + H Ii: (8) 
5 



tin — \ (9) 

for the adiabatic and thermal decoherence regimes respectively, and 

H Ic = -^V(r)[C x a x + C y a y ] } (10) 

where 

C x = c x cos(u t) + c y sm(u t) (11) 
C y = — c x sin(u; t) + c y cos(u t). (12) 

These interaction-picture Hamiltonians have been calculated by using the rotating wave 
approximation, and assuming zero detuning, 8 = ujyi — uj = 0. 
The evolution equation for p admits the formal solution 

p(s, e, r) = Ui(s, e, r)p(s, e, 0)U}(s, e, r) (13) 

where the evolution operator satisfies the equation 

dUi{s,e,r) % 

dr =-fi g /*Mw) (14) 

given the initial condition Ui(s,e,0) = 1. The formal solution of Eq. (jl4|) is 16], 

Ufa^r) = Tiexpi—^dr'Hrir')}], (15) 

where T[ ] represents the time-ordering operator. The state of the system at r = is 
assumed to be described by the density matrix: 

p(s,e,0) = p(s,0)®p(e,0) (16) 

where the environment is in thermal equilibrium, 

p(e,0)=Y[p q (e,0) (17) 
<? 

p q (e,0) = [l-e-^J-^e-^K)^! (18) 
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and 

2 

p(a,0) = X)py(0)P«. (19) 

Since the various terms in the interaction Hamiltonian Hi do not commute, the complete 
evolution operator, Ui(s,e,r), cannot be calculated exactly, except in very special cases. 
By writing the interaction Hamiltonian as a sum of two non-commuting operators, Hj = 
Hi + H 2 , we can expand the evolution operator into an infinite product of exponentials, 

Uj(s, e, r) = e Hf/o T x e -s/o T «^(*) x e -(f) 2 K<ttf* dt'^t)^')] x _ _ _ ( 20 ) 



using the general Baker-Hausdorff theorem 



17| . Assuming that the incremental effects of 



the control and interaction Hamiltonians are relatively small, we approximate the evolution 
operator to the first order in the magnitude of the control pulses, V(r), and the coupling 
strength parameter, e, of the system-environment interaction. 

We consider the evolution operator, Ui(s,e,r), and the density matrix elements for the 
adiabatic and thermal decoherence regimes, separately. The latter case describes population 
changes in the two-level system, and encompasses the evolution of a spontaneously-emitting 
atom, in the quantum feedback scheme proposed by Wang et al. jlfl| . 

A. Adiabatic Decoherence Regime 

The adiabatic decoherence regime describes phase decay of the two-level system. For 
the sake of simplicity, we assume that the control Hamiltonian is proportional to a x , i.e., 

(C x = l,C y = 0), 

Hi e =~V(r)a x , (21) 

and we identify the control Hamiltonian, Hj c , with H%, and the system-environment inter- 
action Hamiltonian, Ha, with H 2 . Therefore, by neglecting the commutator terms in Eq. 
(|27?j) the evolution operator can be approximated by 

U I (s,e,r)^exp{~ f dtH Ic (t)} x ex V {- % - [ dtH R (t)}. (22) 
n Jo riJo 

This approximation is correct to first order in the magnitude of the control and interaction 
Hamiltonians. 

Performing the time integrals in the exponents, we obtain 

t/r(s,e,r) = exp{za x J(r)} x exp{-2a 2 (5_(r)} (23) 



where 



and 



I{r)= l -j^dr'V{r') (24) 



Q-(r) = lJ2( M A- M : a < 1 ) (25) 
q 



M q (r) = ^(l-e^ T ). (26) 



where 



We obtain the following approximate expressions for the evolution operator and its ad- 
joint: 

Ui(s, e, r) = E + iE x o x + E y o y + E z o z (27) 

and 

Uj(s, e, r) = E - iE x a x - E y a y - E z a z (28) 

E = cosh(2g_(r))cos(/(r)) (29) 

E x = cosh(2g_(r))sin(/(r)) (30) 

E y = -sinh(2g_(r))sin(7(r)) (31) 

£ z = -sinh(2g_(r))cos(7(r)). (32) 

The density matrix elements of the reduced two-level system can then be solved explicitly 
by tracing out the environmental modes, 

2 

Pij (s,r) = Tr e {Y,(i\UiPkip(e,Q)Uy l \ 3 )p kl m (33) 

k,l=l 

for i,j = 1, 2, assuming that the environment is initially in a thermal state. We obtain: 

Pn = pn(0)cos 2 / + p 22 (0)sin 2 /-i[p 12 (0) -p 2 i(0)]e^ d cos/sinJ, (34) 

P22 = P2 2 (0)cos 2 / + p 11 (0)sin 2 / + i[p 12 (0) -p 21 (0)]e^cos/sin/, (35) 

P12 = pn(0)e~ 9ad cos 2 / + p 21 (0)e^ d sin 2 / + i(p 22 (0) - pn(0)) cos/sin I, (36) 

P21 = P 2i(0)e- 9ad cos 2 J + p 12 (0)e-^ d sin 2 /- i(p 22 (0) -pn(0)) cos/sin J, (37) 
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where I = I(t), and 



oo 



</.,./ :=</,„/( r) = - / dujG{<jj)(l - cosojt) coth ^ (38) 







is the adiabatic decoherence function obtained by Palma jig ]. The dimensionless constant 
7 depends on the dipole moment of the two-level system and on the Rabi frequency. 



B. Thermal Decoherence Regime 

The thermal decoherence regime describes population changes in the two-level system, 
and encompasses the evolution of a spontaneously-emitting atom. We use the general form 
of the control Hamiltonian in Eq. (JlU)) . and we decompose the interaction Hamiltonian Hj 
into two operators, Hi and H2, which are proportional to a x and a y , respectively. 

By neglecting the commutator terms in Eq. ()20|) . the evolution operator can be approx- 
imated by 

U I (s,e,r)^exp{~J^ dt(H ICx (t) + H Iix (t))} x ex V {- % - J dt(H ICy (t) + H Iiy (t))}, (39) 

where Hj Cx , Hn x and Hi Cy , Hn y correspond to the terms in the control Hamiltonian Hi c and 
the system-environment interaction Hamiltonian Ha, which are proportional to a x and a y , 
respectively. This approximation is correct to first order in the magnitude of the control 
and interaction Hamiltonians. 

Performing the time integrals in the exponents, we get 

C//(s, e, r) = exp{iC x a x I(T) - o x Q_(t)} x exj){iC y a y I(r) + ia y Q + (r)} (40) 

where 

q+(t) = ^J2( m A+ m >«)' ( 41 ) 

MJt) = ^ (1 - e^ 12 -^ T ), (42) 

and J(r), Q_(r) are defined by Eqs. ip^ - lffSj) . 

We obtain the following approximate expressions for the evolution operator and its in- 
verse: 

U T (s, e, t) = E + iE x a x + E y o y + E z a z (43) 
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and 

U}(s, e, t) = E — iE x a x - %E y a y + iE z a z (44) 

where 

E = cos(CJ(r) + zQ_(r)) cos(C v /(r) + Q+(r)) (45) 

£ x = sin(C x /(r) + iQ-(r)) cos(C y I(r) + Q+(r)) (46) 

£ y = cos(a/(r) + zQ-(r)) sin(Cj,/(r) + Q+(r)) (47) 

£ 2 = sm(CJ(r) + iQ_(r)) sin(CV/(r) + Q + (r)). (48) 

To calculate Uj(s,e,r), we have used the fact that Q-(t) and Q+(r) commute in the first 
order approximation of the coupling strength parameter, e. 

The environment modes are then traced out to obtain the elements of the reduced density 
matrix, in the same way as for the adiabatic case (see Eq.([33|)). assuming as before, that 
the environment is initially in a thermal state: 

Pn = i + ^(pn(0)-p 22 (0))e- 2 ^cos(2a/)cos(2C ?/ /) 

-Re{p 12 (0)}e- 29th cos(2C x I) sm{2C y I) - ilm{p 12 (0)}e- 9th sm{2C x I) (49) 
P12 = [i2e{pi 2 (0)}cos(2C„/)+/m{pi 2 (0)}cos(2C x /)]e- ff « h 

+iRe{p 12 (0)} sin(2CU) sin(2C y J)e~ 2sth 

+\(Pn(0) ~ p 22 (0))[e-^ sin(2<7 tf /) - ie" 2 ** sm(2C x I) cos(2C y I)} (50) 
P2i = [Re{p 21 (0)} cos(2C y I) +Im{p 21 (0)} cos(2C x I)]e' 9th 
-iRe{p 21 (0)} sm(2CJ) sin(2C y I)e- 29th 

+ ^(Pn(0) - p 22 (0))[e-^ sm(2C y I) + ie~ 2 ^ sin(2CJ) cos(2C v J)] (51) 

P22 = ^-^(pn(0)-p 22 (0))e- 2 ^cos(2a/) 008(2^/) 

+ J Re{p 12 (0)}e~ 29th cos(2C x I) sia(2C y I) + ilm{p 12 (0)}e- 9th sm(2C x I) (52) 

where the thermal decoherence function, g t h, is given by 

9th 9th(r) = 7 I duj- — ms ^ UJl2 u 3 coth((3 uj/2) exp(-u/uj c ). (53) 

Jo (wi2 - wy 

III. CONTROL STRATEGY 

The control strategy in Ref. Q| was to continually adjust the control pulses in a cus- 
tomized fashion to counter the effects of the decoherence and maintain the density matrix 
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elements unchanged. In this paper, the goal is to drive the qubit to the final target state 
while, at the same time, maintaining quantum coherence throughout the process. The re- 
duced density matrix of the qubit depends on the known decoherence function and the 
applied control, which is determined by equating each element of the actual reduced density 
matrix, pij, in turn, with the corresponding element of the density matrix for the target 
state. 

Since in general, a large jump in the state of the qubit cannot be achieved in one step,the 
qubit is driven between the initial and target states via a number of intermediate states. 
This ensures controllability i.e., the existence of feasible solutions for the proposed control. 

Each control cycle contains 8 control steps, since there are 8 real equations to solve, for the 
reduced density matrix of the two-level system. The sequence of eight real transcendental 
equations is solved in a prescribed order. While this particular order is not essential, we 
shall use it consistenly here, as shown below, to make it easier to follow our strategy. The 
algorithm is applied identically for either the adiabatic or the thermal case, thus we describe 
it for a generic decoherence function, denoted g. 

The intermediate states spanning the trajectory between the initial and target values can 
be chosen explicitly, the simplest way being to interpolate linearly between the initial and 
target states, or alternatively, they can be selected by minimizing the absolute difference 
between the initial and target values for a restricted range of control pulses. 

In the first cycle of eight time steps, we start by driving the first component Phr(O) to a 
new intermediate value, Cur? using the control pulse 1(1), which is a solution of the equation 
Pilr(1) = Piii?(fi , (l), -f(l)) = Cur- After the first time step, all the other components have 
been modified by the effect of the decoherence, g(l), and the effect of the control pulse, 1(1), 
so the values of the matrix elements are given by: 
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Piijs(I) 


= pn«(^(i),/(i))=c 1 ( ;i ? 


(54) 


Pllj(l) 


= p lu {g{l)J{\)) 


(55) 


Pl2i?(l) 


= p 12R (g(l),I(l)) 


(56) 




= pv, t (q(1),I(1)) 


(57) 


P21i?(l) 


= p 21R (g(l),I(l)) 


(58) 


021/(1) 


= p 21I (g(l)J(l)) 


(59) 


P22i?(l) 


= p 22R (g(l),I(l)) 


(60) 


P22/(l) 


= p 22I (g(l),I(l)) 


(61) 



During the next seven time steps, the control strategy is implemented in a similar fashion. 
For example, at the second time step, we consider the equation for pn/(2) = p UI (g(2), 1(1) + 
1(2)). For relatively short time steps, and by the semigroup property, this is also equal to 
([^j(g(l), 1(2)). To determine the control pulse required to set pn/(2) equal to a new value 
we find the control 1(2) which solves the equation pn/(2) = ([^j. After applying this 
control, the values of the matrix elements at the second time step are given by: 

Pij R / i(2)=p mi tg(2),I(l) + I(2)). (62) 

By the end of the first cycle of eight time steps, the state has shifted slightly towards 
the first intermediate state, The final values of the eight components at the end of the 
first cycle become the initial values for a second cycle of eight time steps, which then aims 
to drive the qubit to the second intermediate state, (^ 2 \ using the same method. 

This procedure is repeated until the final target state is reached, and is then maintained 
for as long as required. The total transition time to drive the qubit from the initial to the 
target state depends on the number of intermediate states and the frequency of the control 
pulses. After the initial transient, the controls stabilize and the whole control cycle repeats 
periodically. 
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IV. QUANTUM FEEDBACK CONTROL 



In this section, we shall briefly review the quantum feedback scheme proposed by Wang et 
al. [it]]. The theory of quantum feedback has been developed in the last ten years [3], and 
recently there has been growing interest in using it to cancel the effects of decoherence. This 
can be achieved by using the continuous measurement record resulting from the coupling of 
a system to its environment, to cancel out the effect of the interaction with the environment. 

The system considered by Wang et al. ^(J is a two-level atom, which is driven by 
a resonant classical field of strength a G (—00,00). The lowering operator is denoted 
by er = |1)(2|; its adjoint = |2)(1| is the raising operator, and performs the inverse 
operation. The environment is represented by electromagnetic modes into which the atom 
emits spontaneously at a decay rate 7. The reduced master equation describing this system, 
with no feedback, is derived from the interaction Hamiltonian 

Hi = haa y + h{a^^2g*a q + a^2g q a\}, (63) 
q q 

which is identical to the interaction Hamiltonian considered in Section 2.2 for the thermal 
decoherence regime. Phase decay of the two-level atom, due to atomic collisions or other 
processes, has been ignored in the quantum feedback scheme, i.e. there is no o^-coupling for 
the system-environment interaction. The reduced master equation is given by 

P = -ia[(T v , p] + V[^a]p, (64) 

where the Lindblad superoperator T> is 

V[A]B = ABA 1 - {A^A, B} /2. (65) 

We note that the same reduced master equation can be obtained for the spin-boson model 
we presented in Section 2, in the thermal decoherence regime, using a control Hamiltonian 
proportional to o y . 

The effect of the classical driving field a is to rotate the state of the atom around the y- 
axis, in Bloch space. The Bloch space coordinates are defined by the vector (x, y, z), related 
to the state density matrix by: 

p = ^(I + xa x + ya y + zcr z ). (66) 
13 



The distance from the centre of the Bloch sphere, r = x 2 + y 2 + z 2 , measures the purity 
of the state, with r = 1 for a pure state, and r = for a maximally-mixed state. The 
stationary solutions for the driven two-level atom, 

Xss = 2 i g 2 67 

■y 2 + 8cr 

y ss = (68) 
-7 2 

Zss = 2 i g 2 ' ( 69 ) 

■y 2 + 8a 2 

are limited to the lower half of the Bloch sphere [z < 0), in the x — z plane, and in general 
have low purity, as shown in Fig. 2 of Ref. [l0|. The ground state with no driving is 
the only stationary solution that is a pure state. We note that in general, the interaction 
Hamiltonian due to the classical driving field can be taken as a linear superposition of a x 
and a y , as we have done for our control Hamiltonian Hj c . Experimentally this superposition 
of polarizations can be adjusted by controlling the phase of the classical driving field. 

In their scheme, Wang et al. assume that all of the fluorescence from the atom is collected 
and „ 8 ed in a ho m od y „e detection setup (see Fig. 1 in Re, Q), to — the inte^ence 
between the system and a local oscillator. In the ideal limit of large local oscillator amplitude, 
the point process of photocounts is changed into a continuous photocurrent with white noise. 
The interference is measured by the (normalized) difference between the recorded mean 
photocurrents, 

AI(t) = ^(a x ) c (t) + C(t), (70) 

given that the phase of the local oscillator is set to zero. The subscript c means conditioned, 
as system averages are conditioned on the previous photocurrent record. The last term, ((t), 
represents Gaussian white noise so that, 

C(t)dt = dW{t), (71) 

an infinitesimal Wiener increment defined by 

{dW(t)} 2 = dt, (72) 
E[dW(t)] = 0. (73) 

The conditioning process can be made explicit by calculating the changes of the system 
state in response to the measured photocurrent, using a stochastic Schrodinger equation, 

d\A(t)) =A c (t)\Mt))dt + B c (t)\Mt))dW(t), (74) 
14 



which describes possible quantum trajectories, assuming that the initial state is pure at 
some point in time. The operator A c (t) for the drift term, and the operator B c (t) for the 
diffusion term are both conditioned, as they depend on the system average 

(<r x )c(t) = (Mt)\**\Mt))- ( 75 ) 

Of course, from the point of view of ensemble averages, the two-level atom still obeys the 
same master equation given above in Eq. (JMJ). 

The difference between the recorded mean photocurrents, A/(t), is used to control the 
dynamics of the state in the x — z plane. The atom-environment dynamics is altered by 
adding a photocurrent feedback, XAI(t), to the amplitude of the classical driving field, using 
an electro-optic amplitude modulator. The effect of this feedback term is to change the drift, 
A c (t), and the diffusion, B c (t), operators, such that the feedback-modified master equation 
describing the quantum feedback scheme becomes: 

p = -ia[a y ,p] +V[ v ^a - i\a y ]p. (76) 

The two-level atom can be driven to, and subsequently maintained in an arbitrary pure 
target state \6), in the x — z plane, 

|0}=cos(^|2)+sin(f)|l), (77) 



x 2 / V 2 , 

by ensuring that only the phase of \9) changes under this evolution: 

\6) + d\6) = e ldx \6). (78) 

In terms of the stochastic Schrodinger equation, this implies that the deterministic (drift) 
and noise (diffusion) terms are separately proportional to \8), 

A c (t)dt\e) oc |0) (79) 
B c (t)dW(t)}\6) oc \9), (80) 

which can be achieved by setting the values of the driving and feedback parameters to: 

a = 7/4 sin 9 cos 8, (81) 
A = -^5(l + cos0). (82) 
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A stability analysis of the scheme shows that the pure state \9) can be achieved for all 
9, with the exception of 6 = ±7r/2, which correspond to equal superposition states on the 
equator of the Bloch sphere. The driving and feedback parameters for these special values 
are degenerate, thus the state is driven to a mixture of |~) and | — ~). 

The quantum feedback scheme described above relies on two major assumptions, that are 
at present impossible to achieve experimentally, namely: (i) perfect detection of spontaneous 
emission, and (ii) zero feedback delay time. 

As the efficiency of the homodyne detection of resonance fluorescence from the two-level 
atom is decreased from perfect efficiency, the performance of the scheme degrades very 
quickly, except for target states that are close to the ground state. In general, detection 
with non-unit efficiency gives a stationary solution that represents a mixed state, rather 
than the desired pure target state. Wang et al. have shown that the deviation from a pure 
target state is noticeable even for relatively high (95%) efficiencies, especially for states close 
to the equal superposition states at the equator of the Bloch sphere, which are inherently 
unstable. 

Similarly, in order to successfully achieve a pure target state, the feedback delay time 
t has to be much shorter than the time scale corresponding to the spontaneous emission 
rate. This limit corresponds to the Markovian regime, which allows the evolution equation 
to be written in the form of a master equation. If the feedback delay time becomes too 
large, r > 0.027 -1 , the performance of the quantum feedback scheme suffers a degradation 
which is qualitatively similar to that observed for inefficient detection. For example, if we 
try to stabilize the atom in the excited state, in the non-Markovian regime, the purity of 
the average state is given by r w 1 — A'-fT. Of course, the deviation from a pure target state 
is much more noticeable for states close to the equal superposition states at the equator of 
the Bloch sphere. 

We have verified that the quantum feedback scheme works as intended, in the regime of 
unit detection efficiency and zero feedback delay, by simulating typical quantum trajectories 
for the feedback master equation in Eq. (J7HJ). The transition time, t , required to drive the 
atom from an initial state to the target state, in the x — z plane, is roughly proportional 
to the distance between these states, as measured by the difference in the angle 9. The 
longest single-transition times (~ (4 ± 1)7 -1 ) correspond to driving the atom between two 
states that are diametrically opposite on the Bloch sphere, e.g., from the ground state to 
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FIG. 1: The two-level atom is driven from the ground state, |1) = \ir), to the excited state, \2) = |0). 
The single quantum trajectories (crossed points) shown in (a) and (b) represent the two different 
types of paths that the state of the atom may follow. The two paths correspond to opposite 
directions around the Bloch sphere, passing through the (a) symmetric, or (b) anti-symmetric 
superposition of the ground and excited states. The fidelity of \8) compared to the target state |0) 
is also shown (solid line) on the graphs. The transition time, to, defined here as the time taken to 
reach a state fidelity |(0|6)| 2 > 0.99 is (a) t = 4.07" 1 and (b) t = 3.67 -1 . 
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FIG. 2: The non-zero density matrix elements of a qubit driven from the ground state, |1), to 
the symmetric superposition state, (|1) + i\2))/y/2. All of the density matrix elements show a 
monotonic transition to their target values. This driving used 6 intermediate states, with only 
one control cycle allocated for each intermediate state, and the maximum allowed strength of each 
control pulse was set by I ma x = 0.1. Note that the elements change in a noticeably discrete way 
during the transition period, and the fidelity of the final state is greater than 0.99. 



the excited state, as shown in Fig. [TJ 



V. RESULTS AND COMPARISON 



We first summarize the results of our open-loop control scheme, and then, whenever 
warranted, present a comparison between our results and the results of the quantum feedback 
scheme proposed by Wang et al. for a two-level atom. 
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A. Results for Open-loop Control 

First, we present the open-loop results for the evolution in the adiabatic decoherence 
regime, which corresponds to phase decay of the two-level system. Using a control Hamil- 
tonian proportional to a x , the qubit can be driven to any target state on the surface of the 
Bloch sphere, in the y — z plane. Only a small number of intermediate states are required 
to obtain robust control. The graphs in Fig. 2 illustrate the driving of the qubit from the 
ground state |1) to the symmetric superposition state (|1) + i\2))/\/2. The qubit was then 
driven to the excited state, and back to the ground state via the conjugate anti-symmetric 
superposition (|1) — i\2))/y/2, to complete a continuous circle around the Bloch sphere. 

It is also possible to jump directly from the initial to the target state in a transition time 
as small as one control cycle, if the control pulses are strong enough. However, this can 
break the assumption that the magnitude of the control Hamiltonian is small, and moreover 
gradual driving improves the smoothness of the transition, albeit at the expense of longer 
transition times. 

We have not simulated the control for the adiabatic decoherence regime using other 
control Hamiltonians, as there are no analogous results for the quantum feedback control of 
phase decay, either by itself, or in conjunction with spontaneous emission. 

Next, we present the open- loop results for the evolution in the thermal decoherence 
regime, when the control Hamiltonian is proportional to a y . This corresponds to the en- 
semble evolution of the spontaneously-emitting two-level atom (without feedback) in the 
quantum feedback scheme. The numerical simulations used to implement our control strat- 
egy from Section 3, have assumed the ideal scenario: the decoherence function is known 
exactly, and the time step between control pulses for the open-loop control is small relative 
to the decoherence rate. 

The graphs in Fig. 3 illustrate the driving of the qubit from the ground state 1 1) to the 
excited state |2) via the symmetric superposition (|1) + \2))/y/2, cf., Fig. 1(a) for the same 
driving in the quantum feedback scheme. One hundred intermediate states were used in this 
simulation to give a very smooth transition in a time t ~ 1250At, where At is the time 
interval between control pulses. The qubit was then driven back to the ground state via 
the anti-symmetric superposition (|1) — |2))/y / 2, to complete a full circle around the Bloch 
sphere. We note that returning to the ground state with exactly the same phase, actually 
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FIG. 3: The non-zero density matrix elements of a qubit driven from the ground state, |1), to the 
excited state, |2). The population elements, Re{pu} and Re{p22}, show a monotonic transition, 
whereas the coherences, Re{p\2} and Re{p2i}, increase to | as they approach the superposition 
state on the equator of the Bloch sphere, and then decrease back to zero. This driving used 100 
intermediate states, with only two control cycles allocated for each intermediate state, and the 
maximum allowed strength of each control pulse was set by I ma x = 0.01. All of the density matrix 
elements change smoothly in the transition period, and the fidelity of the final state is greater than 
0.99. 



requires two complete circles around the Bloch sphere. 

In general, our simulations of open-loop control in the thermal decoherence regime, show 
that the qubit can be driven to any pure state in the x — z plane of the Bloch sphere, 
using a control Hamiltonian proportional to a y . Inevitably, the fidelity of the target state 
thus achieved fluctuates with each control cycle, but it stays close to one, given that the 
maximum allowed strength of the control pulses is relatively small. 

The smoothness and length of the transition to the target state is determined by the 
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FIG. 4: The non-zero density matrix elements of a qubit driven from the ground state, |1), to the 
excited state, |2). Only 10 intermediate states were used for this driving, with two control cycles 
allocated for each intermediate state, and the maximum allowed strength of each control pulse was 
set by Imax = 0.1. All of the (non-zero) density matrix elements display large fluctuations in the 
transition period, and the fidelity of the final state is less than 0.99. 

number of intermediate states, as well as the maximum allowed strength of the control 
pulses. The graphs in Fig. 4 illustrate explicitly the same driving of the qubit as in Fig. 3, 
but using only ten intermediate states, and allowing the control pulses to be up to ten times 
stronger than before. We note that in the thermal decoherence regime, there is an absolute 
minimum number of intermediate states, of the order of ten states, to ensure successful 
control of the qubit, for the longest type of single-transition. 

The transition time from the initial state to the target state depends not only on the 
number of intermediate states, and the maximum allowable strength of the control pulses, 
but also on the magnitude of the time step, At, between control pulses. The minimum times 
we observed for the longest single-transition of the qubit e.g., from the ground state to the 
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excited state, are of the order of ten control cycles, for the type of low-quality results shown 
in Fig. 4. Of course, the transition times grow longer as the performance of the open-loop 
scheme improves, e.g., the driving shown in Fig. 3 requires a transition time of the order of 
100-200 control cycles. 

We noted earlier, that our control Hamiltonian, and the equivalent interaction Hamilto- 
nian due to the classical driving field in the quantum feedback scheme, can be taken as a 
linear superposition of a x and cr y , see Eq. fllOj) . In general, if the control Hamiltonian is 
written in terms of the angle 0, 

H Ic ((fi) oc (cr x sin (ft + a y cos (f>), (83) 

the stationary solutions of the evolution, 

4cry cos (j> 



7 2 + 8a 2 



(84) 



-4«7sin v 

V" = 2 i g 2 85 

(86) 



ss 7 2 + 8a 2 ' 

are contained in a plane S^, and lie on a curve parametrized by the ratio, a/7, of the driving 
field and the decoherence rate. 

The plane S^, defined by xsmcf) + ycoscj) = 0, always contains the z-axis (i.e., always 
includes both the ground and the excited states of the atom), and is rotated by an angle 
around the z-axis in Bloch space, with respect to the reference x — z plane (S^ = o). Thus the 
effect of the driving field corresponding to H Ic (<fi) is to rotate the state of the atom around 
an axis which is normal to the plane S$. 

It is evident that if we employ different control Hamiltonians, Hi c ((j)), the stationary 
solutions for the evolution of the spontaneously-emitting two-level atom (without feedback) 
are not restricted to the x — z plane. This result also applies to the quantum feedback scheme 
of Wang et al. and our open-loop control in the adiabatic case, enabling us to stabilize a 
qubit in any pure state (except for equal superpositions) on the surface of the Bloch sphere. 

In general, the atom can be driven reversibly to any arbitrary pure target state on the 
Bloch sphere, by a combination of two control Hamiltonians. The first Hamiltonian rotates 
the initial state to the ground state, and the second Hamiltonian rotates the ground state 
to the desired target state. To return to the initial state, the order and the sign of the 
Hamiltonians need to be reversed. 
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FIG. 5: The non-zero density matrix elements of a qubit driven from the ground state, |1), to the 
excited state, |2). The population elements, Re{pu} and Re{p22}, show a monotonic transition, 
whereas the coherences, Im{p\2} and Im{p2i}, increase to ^ as they approach the superposition 
state on the equator of the Bloch sphere, and then decrease back to zero. This driving used 100 
intermediate states, with only two control cycles allocated for each intermediate state, and the 
maximum allowed strength of each control pulse was set by I ma x = 0.01. All of the density matrix 
elements change smoothly in the transition period, and the fidelity of the final state is greater than 
0.99. 



To demonstrate this idea, we have also simulated driving the qubit in the thermal deco- 
herence regime, using a control Hamiltonian proportional to a x . The same principle applies 
if the control Hamiltonian is a linear superposition of a x and a y , though the computation 
of the evolution operator Uj and the density matrix elements is more intensive. 

The graphs in Fig. 5 illustrate explicitly the driving of the qubit from the ground state 
|1) to the excited state |2), using a control Hamiltonian H' Ic oc a x , c.f., Fig. 3 for the 
same driving with H Ic oc o y . The state of the qubit is confined to the y — z plane of the 
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Bloch sphere (and passes through the symmetric superposition state + i|2))/v2), thus 
the effect of the driving field is to rotate the state of the atom around the x-axis. One 
hundred intermediate states were used in this simulation to give a very smooth transition 
in a time t Q « 1250At, where At is the time interval between control pulses. The qubit 
was then driven back to the ground state via the conjugate anti-symmetric superposition 
(|1) — i\2))/y/2, to complete a continuous circle around the Bloch sphere. 

To illustrate how a combination of two control Hamiltonians can be used to achieve any 
arbitrary pure target state, suppose that we are given the initial state if)i n iuai — (|1) + |2))/ v2, 
in the x — z plane, as shown in Fig. If we are asked to drive the qubit to the target 
state iptarget = + i|2))/ a/2, in the y — z plane, our control strategy is as follows: first 
use the control Hamiltonian Hj c oc a y to drive the qubit to the ground state, along the 
curve determined by the intersection of the Bloch sphere with the x — z plane, and then use 
the control Hamiltonian H' Ic oc a x to drive the qubit to the target state, along the curve 
determined by the intersection of the Bloch sphere with the y — z plane. 

B. Comparison between the Open- loop and Quantum Feedback Controls 

The open-loop evolution in the thermal decoherence regime, is identical to the ensemble 
evolution of the spontaneously-emitting two-level atom (without feedback) modeled by Wang 
et al. when the control Hamiltonian is proportional to a y . After quantum feedback is 
introduced, the difference in the dynamics between the two schemes becomes evident: the 
feedback directly changes the effective interaction of the atom with its environment, whereas 
in the open-loop scheme, the control Hamiltonian, which is equivalent to the classical driving 
field a in Eq. (|63|h affects only the atom. 

Both the open-loop and the quantum feedback control schemes presented in this paper 
aim to drive a qubit to an arbitrary target state, and then maintain that state. For simplicity, 
we have restricted our attention to pure states, though in principle the initial and target 
states may be mixed states in both schemes. The general case will be investigated in future 
work to check our conjecture that it is only possible to drive to a target state that has a 
purity less than or equal the purity of the initial state. 

The initial conditions are not important for the quantum feedback scheme, whereas with 
open-loop control we must assume prior knowledge of the initial state. This requirement 
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FIG. 6: The effect of two orthogonal control Hamiltonians for pure states on the surface of the 
Bloch sphere. The qubit is driven from the initial state ipinitial to the ground state |1) (solid arrow), 
using a control Hamiltonian Hj c oc a y . Next, a control Hamiltonian Hj c oc a x is used to drive the 
qubit from the ground state to the target state ^target (broken arrow). 

can be relaxed for our scheme, if we are able to reset all initial states to the ground state, 
so that we have a fixed starting point for the qubit. Finally, neither scheme can be used 
to maintain an unknown state such as that obtained from a quantum computation. Both 
schemes can only drive the atom to a known target state. 

If the evolution in our model is the same as the ensemble evolution in the quantum 
feedback scheme, the open-loop control can drive the qubit between any two pure states 
in the x — z plane (see the results illustrated in Fig. 3). Therefore, our simulations have 
demonstrated that open-loop control can effectively achieve the same result as the quantum 
feedback scheme. 

We discuss in turn the (i) efficiency, (ii) feasibility, and (iii) robustness of the open-loop 
and quantum feedback schemes: 



25 



1. Efficiency. 



The final state actually achieved by the open-loop control scheme has a fidelity which 
fluctuates during each control cycle, but remains very close to one with respect to the target 
state, given that the maximum allowed strength of the control pulses is relatively small. The 
transition times in our scheme are much shorter than the equivalent times for the quantum 
feedback scheme, if the rate of the control pulses is higher than w IOOO7. 

2. Feasibility. 

The quantum feedback scheme relies on the assumptions of perfect detection and zero 
delay time. In the case of inefficient detection, or in the non-Markovian regime (non-zero 

n 

delay), the final state is a mixed state [10]. The deviation from a pure target state is 
noticeable even for relatively small violations e.g., 95%-efficiency or Tp > 0.027 -1 , especially 
for states close to the equal superposition states at the equator of the Bloch sphere, which 
are inherently unstable. In terms of practical implementation, it is currently impossible to 
satisfy these assumptions due to technological challenges. 

On the other hand, the requirements for the open-loop control scheme are in principle less 
stringent. The time step between control pulses needs to be small relative to the decoherence 
rate, which can be satisfied in a practical implementation, by choosing an atomic transition 
with a decoherence rate that is low compared to available control pulse rates. The strength 
of the control pulses and the strength of the system-environment interaction also need to be 
small relative to the free Hamiltonian, in order to validate the first-order approximation of 
the evolution operator obtained in Eqs.(j2HJ) and PU|) . The other major assumption is that 
the decoherence function is known accurately, which is essentially true for a single two-level 
system such as a trapped atom. 

3. Robustness. 

The stochastic model used in the quantum feedback scheme describes the evoulution of 
single quantum systems using quantum trajectories with an extrinsic noise source. On the 
other hand, the spin-boson model assumes knowledge of the average decoherence function, 
and a deterministic evolution, which corresponds to an ensemble of quantum systems. 
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The stochastic model is a more realistic portrayal of the quantum dynamics of a single 
spontaneously-emitting two-level atom than the generic ensemble model we have employed. 
Thus, in principle, the quantum feedback scheme is always going to perform more robustly 
than any open-loop scheme, because it is inherently able to use feedback to control a single 
qubit. Following the transition period, the effect of random fluctuations on the state of a 
single qubit are effectively canceled by the feedback. 

The spin-boson model uses an average decoherence function describing the behavior of 
an ensemble of qubits. Our scheme cannot respond adaptively to random fluctuations for 
a single qubit, which will result in much greater deviations from the target state, and the 
open-loop control may fail completely, if the fluctuations are too large. 

In conclusion, we presented an open-loop control scheme that drives a qubit to an arbi- 
trary target state with high fidelity, while maintaining quantum coherence throughout the 
process. 

We have also shown that by using different control Hamiltonians, which are proportional 
to a linear combination of a x and a y , the quantum feedback scheme of Wang et al. can 
be extended to stabilize the qubit in any known pure state, except for equal superposition 
states. 

The performance of our scheme compares favorably with the performance of the quantum 
feedback scheme proposed by Wang et al, with respect to target fidelity, and transition time. 
In principle, both schemes can drive a qubit to an arbitrary pure target state, with high 
fidelity, in a time of the order of lOy -1 or less. 

From a practical viewpoint though, while the quantum feedback scheme has the capability 
to react to sudden fluctuations for a single qubit, its performance suffers from the strict 
technical requirements for zero feedback delay and unit detection efficiency The open-loop 
control scheme has the additional advantage that the transition time can be lowered by 
simply increasing the rate of the control pulses. 
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